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ABSTRACT 


This paper presents the development of a second order theory for trajec- 
tories in the vicinity of the lunar libration point L ? . This development is 
based on a four-body model including the sun, earth, moon, and a satellite, 
all assumed to move in the same plane. 

As a result we obtain a system of two simultaneous second order differ- 
ential equations with time dependent coefficients. Some selected solutions 
of interest are derived here, including a first order periodic solution and a 
second order quasi-periodic solution. Various trajectories around L are 
plotted and computations of the velocity, acceleration, range, range rate and 
flight path angle are presented. 
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1 . 


INTRODUCTION 


The general three body problem is known to admit only one solution, the 
one found by Lagrange. This solution is satisfied at the five libration points 
of the earth-moon system. The libration point, L ? , located at the far side of 
the moon is of special interest to scientists lately because it could provide 
an "anchor" for a communications satellite behind the moon. 

However, the three body model yields only a rough first order approxima- 
tion to the motion of a satellite in the vicinity of L s . Even this solution 
demonstrates that is an unstable "anchor" or, in other words, a satellite 
will not stay in orbit around L 2 unless forced to by an onboard variable 
thrust engine. 

It is precisely for this reason that we are interested in finding out an 
improved approximation to the motion about 1^ . The relevant factors that do 
not appear in this earth-moon model include the gravitational field of the 
sun, the oblateness of the earth, the eccentricity of the moon 1 s orbit, and 
the inclination of the moon T s orbit to the earth's equatorial plane. Two other 
factors that are also excluded are the pressure of solar radiation and meteo- 
roid disturbances. Of all of these external perturbative forces, the gravita- 
tional field of the sun is the most important. 

In this analysis we will consider first and second order effects of the 
sun on the motion of a body around L 2 ; these effects introduce a non-homogeneous 
forcing function. We construct a four body model consisting of the sun, earth, 
moon and a satellite stationed initially at 1^ or its immediate vicinity. The 
sun and the moon are assumed to move in circular coplanar orbits with respective 
constant angular velocities 0 and co. By assuming uu to be constant, we clearly 
neglect the eccentricity of the moon's orbit. An expression for a generalized 
acceleration, including solar perturbation, is then developed. Its components 
are the force functions of the equations of motion. 

The complete solution is given in terms of first and second order solu- 
tions derived by the method of regular perturbations. Furthermore, the initial 
conditions are chosen in such a way as to eliminate the dominant unstable 
contribution. This effect can be implemented in practice during the injection 
of the satellite into its orbit around . One has to specify an injection 
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location and then the injection velocity. Although the second order solution 
is unbounded in time, its rate of growth is small and can be corrected by an 
onboard engine. Various computations such as velocity, acceleration, pertur- 
bative acceleration, range, range rate, etc., are presented here as related 
to the trajectory. These computations, which include the effect of solar 
perturbation on the motion, give sufficient information to design a mechanical 
correction for the unstable effect of the motion. 

It should be noted that a first order analysis with solar perturbation has 
been done before for a satellite near the collinear libration points, but, as 
our paper will show, second order effects in solar perturbation are of great 
importance since is an unstable point. Thus, these effects must be included 
to give a good analytic approximation to the motion, as was done here. 
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2 . 


DERIVATION OF THE EQUATIONS OF MOTION 


The following derivation is for a 3-body system in 2 dimensions, describing 
a planar motion. The effect of solar perturbation will be added later. 



Let (Xg, Yg) be a set of inertial coordinates referenced at the earth's 
center. In terms of the rotating coordinates (X, Y), centered at the earth- 
moon barycenter, the equations of motion for a satellite, P, due to the gravity 
field of the earth and the moon are: 


X 


2”jY = x 3 X 




(X-Xi) 

V 


(X-Xo) 

- ^ 


(la) 


5 + - » 2 ir - nj ^ ^ 


(lb) 


where ao denotes the rate of rotation of the earth-moon system, and the terms 2cnX, 

2uuY are the Coriolis accelerations, whereas the terms (jo 2 X, cju 3 Y are the centrifugal 
accelerations . 

At the five u libration points 11 the right hand side of equations (la) and (lb) 
is identically zero, and thus the solutions at these points are X = constant, Y = 
constant. Let the coordinates of the Lg libration point be given by X = X c , Y = Y c . 
Then in order to study the small motion near (X^, Y ), let X = X^ + x and Y = Y c + y 
in (la) and (lb), and expand the r.h.s. about (X^ , Y^) . If we expand only up to 
and including the first two terms of the Taylor series, we will obtain the following 
set of linear differential equations in (x, y) , centered at L 2 . 


- 3 - 



x - 2uuy - (1 + 2A) o)®x = 0 


(2a) 


where A = 


y + 2ujx - (1 - A) «J 2 y = 0 

i r 


^ + ^ -j 


n E + ^ L (1 + p) 3 P 3 


L bl 


E i 


M 


Ls 


EM 


P r 


EM 


Figure 2. Location of L 2 


(2b) 

(3) 


To include solar perturbation we must develop an expression for the accelera- 
tion of the satellite P, relative to the acceleration of the libration point L ? . 
That is: 


r = r - r 
PLg PE LgE 

Then the equations of motion will read: 


(4) 


x - 2uuy - 

(1 + 2A) uj 3 x = a' x 

(5a) 

y + 2<yx - 

(1 - A) w 2 y = a y 

(5b) 


where a^, a denote respectively the x and y components of r. 
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So we proceed to develop expressions for r and ? in terms of known parameters. 

r E Lp E 



Figure 3. 4-Body Configuration 

The acceleration with respect to earth, of the moon, M, due to the gravita- 
tional attraction of the earth, E, and the sun, S, is given by: 


~ = - 4 M ^ _^S - ^S 

EM ' " EM ' rEM r Ms rMS " r£S 


( 6 ) 


The only vector in (fc that is not known directly is f 

MS 


But r s f - t 

MS ES r EM 


(7) 


and r MS r ES ‘ 2r ES * r EM + r EM r ES [ 3 


2f • r r 

ES EM , EM 

•d t • 


ES 


EM 1 

r 2 J 

ES 


Using the binomial expansion and retaining only second order terms in , we 

obtain: 


r MS ~ r ES _ 1 + 3 r 5 


r ES’ r EM 3 r EM i 15 ^ r ES* r EM^ 


ES 


2 r ^ 2 

ES z 


ES 


] 


( 8 ) 
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Substituting (7) and (8) into (6) and rearranging terms, retaining only fourth 

order terms in — t- : 
r ES 


EM 


^ ^ - , bs r r , 3 ( V ? em> - 

r EM ^ ' 


ES 


r s 

3 EM. 

2 Ts*™ 


(r • r 

15 ^ ES EM'’ - 1 

"T- —i r — 


ES 


'ES 


3(r ES >r EM ) 

r ls 


EM 


(9) 


Now 


, r L E = (1 + D ) ^ EM as can be seen from Figure 2, 


then r^ E (1 + o) r EM 


( 10 ) 


and : 


u n + \ tin + 

~ ( o) r| M 'EM 


r 

(1 + 0 ) 

ES 


- , 3(? ES* r EM ) * 

r„.. + -3 r 


EM 


ES 


ES 


3(r ES* r EM ) ^ 

? r 

r * 

ES 


( r iTO * r T?Tur) 


3 EM - ,15 V ‘‘ES 'EM'' - 

r__ + — -f r 

ES 


EM 2 r jf ES 2 

ES 


ES 


] 


( 11 ) 


The acceleration of the satellite, P, with respect to the earth, due to the 
attraction of the earth, moon, and sun, is: 






Jk* + _^s , 

r E 3 M r PS r?S ‘ r ES rES 


( 12 ) 


where r is the distance vector from the earth to the satellite. To obtain a 

Jl Ij 

good approximation for r we perturb the satellite along the earth-moon axis by 

L 1-j 

a distance x measured from L p . The motion is not perturbed in the normal direc- 
tion since the dominant gravity effect is along the earth-moon axis. Since the 
motion is confined to the neighborhood of L P , due to the linearization of the 

equations of motion, x is small compared to r ; but the inclusion of x gives rise 

EM 

to a few large terms in the equation of motion as will be seen now. 

Thus let ? pE = (1 + n ) ? EM + xf (13) 


where 


and 


It 

r 

MP 



denotes a unit vector along 

then follows that: 

= or + xi 
EM 

= r - r 
ES PE 


the x-axis. 


(14) 

(15) 
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From (13) , (14) and (15) we now obtain the corresponding scalars: 


.3 

“PE 


r~_ = (1 + or 


r MP 


2 2 
P r 


EM 


PS 


r. 

= r 2 1 + — 

ES L r. 


r 3 

fl + 

_i x_ + 

1 

EM 

L 1 

1+ P r EM 

(l + o ) 2 

. 2 

X 

, 1 x 2 “j 


4- — 

— 

+ t rr i 


p 

r EM 

p r EM 


r 2 
PE 

2(r 

• r ) 

ES PE ; 1 


pr ■ 

’ 

7 s J 


r ES 


ES J 



< 2 

r 2 

EM 


Applying the binomial expansion and neglecting second-order and higher terms in 

x/r , we obtain: 

EM 


..3 

PE ^ 

(1 + 

p)- 3 

r ™[> 

L — I 

1+p r J 

y EM 



.-3 

_3 

r- 3 r 

i - i 

— 1 



MP ~ 

P 

EM L 

P 

r EM - 



, _3 

T-“ 3 

r, 

r 2 

3 PE 

, 3(? es‘ ? pe> , 

15 

^ r ES’ r PE* 1 

PS 

ES 

L 1 

2 r 2 
ES 

'Is 

2 

C/T 

1 


(16) 

(17) 

(18) 


Substituting (18) in (12) and neglecting fifth-order and higher terms in 1/r 


ES 


we obtain for r 


PE 


it. 


PE 


E 

~3~ 

PE 


“MP 


a~ 

EM 


r +■■ r - “ r + — n — I - r 
r„-„ PEC HP C EH r|g L PE 2 r 


3 r PE 3 ' r ES' r pE^ - 

r„„ + 5 r 


ES 


ES 


ES 


3( ' r ES ,r PE- ) - + 15 ^ r ES* r PE^ - 

r^I r PE 2 rC r ES 


ES 


ES 


(19) 


ES 


Substituting (11) and (19) in (4) , we get the expression for the net acceleration 
of the satellite P with respect to L p : 

X - - ^ S. + 4 ^ + 0 + »> *EH ?EM + *ES L- 

- 3 (3 + c, ? pE ♦ 3 <l +0 , 


ES 


ES 


ES 
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( 20 ) 


(r *r ) 2 
15 v ES PE ; - 

• 4 ES 


2 r. 


ES 


T d-o) r 4 


( r ES‘ r EM> _> 1 

J 


ES 


ES 


Let r 


ES 


^S 1 + y ES ^ 


r = r l 
EM EM 


then (tgg't^) ^s 1 ] 


EM 


^ES* r PE^ ~ (1+P ^ ^S^EM + X ES X 

Using these relations in addition to (13), (14), (16) and (17), we obtain a 
simplified version of (20): 


*-X 

r 


PL, 


t 1 ■ & 4 ] [ <1+ ^ 

.3js_-]r = ._*■-> hi 

P r„ w L p 


f - + xI ] + t 1 


'EM 


'EM 


M-E + Mv. 

EM J rjT r EM + (1+o) - 3 r 


r_ + xi 


r EM EM 


p s r f 3 r Es 2 r , . . , , x 1 . , fEs v 

-3- i- XI - - — s- r p (1+p) + 2 (1+n) - I + 3 -5- x^_x 

r ES 1 2 r ES EM L r EM J r ES 


+ pr [-P d-o) ^ s r m r EM - (l+ P )x ES r EM xi + (l+o)* ES *? EM _ 


ES 


+ f [° <1+P) + 2 <1+0> ’'eVem*] } 

Lo 


( 21 ) 


At x = 0 the contribution to ? of the moon and the earth is zero since a 

PL^ 

libration point in the three body system is a point where the sum of the forces is 
zero. That is : 


at x = 0 r. 


PL. 


(1+P) 


In- 


earth = - 


moon 


(l + p) 


r ™ r “ 


hi 

75 


r 7 " r EM 
EM 


% 

r. 


EM 


EM 


^E + ^ 


r 3 
EM 


r EM ° 


( 22 ) 
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With (22) taken into account, (21) now becomes: 


PL C 


_1_ r r r- r EM ~] 

U 3 ' <l+o)» j L ‘ r m J * r| s 


{-xi - | ^ r 2 [p (1+p) + 2 (1+0) ] + 3 

^ z r ES r EM " 


ES 

c 3 

ES 


X ES 


x 


+ r| s L"° (1+0 ^ X ES r EM r EM ' <1+0 ' X ES r EM Xt + (1+P ' ^Es'^EM 


+ T [° (l +> X E S r EM + 2 « 1+3 > X e“s‘eM"J J 


r xl 1 


(23) 


With x = 0 the above expression for r reduces to that developed by 

XT Lp 

F. T. Nicholson (ref, 4 ), However, in our case first and second order terms 

in solar perturbation appear in the acceleration, and the contributions of the 
earth and the moon are also included. 


In accordance with eqs.(5a) and (5b) we are looking for the rectangular 


components of r. 


PL. 


Let 

r PL P = 

= a i + a j 
x y 


( 24 ) 

and 

break 

r into its 

PL p 

components to 

obtain a and a . 

x y 



2 r ^e 

tVj x + r 

-* - | 7P- r * [0 (1+p) + 2 (1 +D ) ] 

2 r ES EM L r EM 


a ~ 

X 

L < 1+ °> 3 

r ls 


, - ES , 

+ 3 73- x + 
ES 


4 -r. 

Is L 


(1+p) x r 


ES EM- 


15 X ES | 


2 r 4 L p (1 EM 

to 


+ 2 < 1+ “> ^sWl) 


Rearranging and neglecting terms in x of order and higher, we obtain the 

r E s 

final expression for a^. 


[7 


M*r 




M 


x ~ rjf, L(l+p) 

£.n 


T - ^ I * * I » < 1+ °> r E '„ ( 5 (^) -3 (^) 




S a 


,*ES s 




ES n 




^ r f is 

"ES " r ES 


+ r# 3 (-^) 2 -1 

^ T?n — ■n n 


(25) 
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Similarly: 


r 

r ES 


3 J ES 3 r ... . 

— r | M L p(1 p ) 


2 r: 


ES 


+ 2 r + [p( 1+ P> x |s r QI + 2 (l+P ^ *ls r EM J 1 

Cj O 


Rearranging as previously: 


1 Pd*o) ^ ^ L 5 < r* ) s < r 5 ) 

2 r Eg EM L r Eg r Eg 


(^) J 

ES 


+ 3 


ES 


r ES 


) ( 


ES 

f ES 


) x 


(26) 


As can be seen, contains only solar terms. This is an expected result 
since we assumed perturbations only along the x-axis. As stated before, the 
earth and the moon do not exert large forces on the satellite normal to their 
axis as long as the satellite is in the neighborhood of 1^. We simplify the 
equations by introducing the following constants. 


K - 2 [_^E _ 

Kl L(i+o> 3 ^ 


'EM 


7 


O M'O 

K 2 = f P (1+p) 7T- r 2 




r E 4 S EM 


M< c 


K-, = 


ES 


where 


M„ 


m, - k 


E *E + \ 


= (ju 2 r 3 — — 

EM Me + m m 


^M 


= K 


\ + m m 

M„ 


= $ 




= K 


S 


em Me + Mm 
= ^ r EM M e 4 M m 


(27) 

(28) 
(29) 


and — = 388.9237 
r EM 


Thus 


„ P r V ( W V ( W 
K i " 2- |_ (1 + p) 3 " p 3 J 


(30a) 
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(30b) 


K a = 


3 

2 


P 


<1 + P) 


^ r E^S /(M E + V 
(388.9237) 4 


f '^ M S /(M E + V 
(388 .9237) 3 


and 


. V<V +M M> , V®E + V 
A “ — (i+op — — 7 


(30c) 

(30d) 


We adopt the units of kilometers and days so that most of our computations will 
be of the order of one. Furthermore, since we are dealing with long missions, 
the units of days are quite appropriate. The constants used are defined below: 


V (m e + V = °- om 


V (M S + V = °* 9879 

m s /(m e + V = 328 ’ 43 ° 


o = 0.167832 

u) = 0.22997 radians /day 

(f) = -0.2128 radians/day 

r__ = 384,752.7 kilometers 

A = 3.17979 

K-l = -0.20512 radians 3 /day 8 

K s = 0.8587319 x 10 -1 (kilometers) /day 3 

K 3 = 0.29525 x 10 ' 3 radians 2 /day 3 


Now, suppose that the sun moves in a circular orbit, coplanar with the 
moon's orbit in the (x, y) plane, with an angular velocity $. Then: 


= r ES C ° S<4t 

(31) 

y ES = r ES Sinrtt 

Substituting (31) into (28) and (29) we obtain expressions for the 
acceleration in terms of x and t. 

a x w KjX + Kg [5 cos 3 0 t - 3 costft] + K 3 [3 cos 3 0 t-l] x 
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a ^ K P 

y p 


—I 

5 cos^t sin^t + sin^t j + 3 K 3 |^cos(/>t sin^t^ x 


A further rearrangement leads to a simpler mathematical form with double 
and triple angles. 


x 


KiX + - K P 


3 costft + 5 cos3^t | + — K 3 1 + 3 cos2<5t j 


(32a) 


1 r ~j 3 r- 

K P j^sin^t + 5 sin3$t | + j K 3 ^sin2$t 


(32b) 


The resultant acceleration is : 


a = a(x,t) = Va 2 + a 
x y 

at x=0, which coincides with the location of L 0 : 


(33) 


Ks 


26 +8cos 2 0t + 30 cos 3 $t cos^t + 10 sin3$t sintftj 1 /* 


(34) 


and the maximum acceleration at x=0 is 


a = 2 K 0 = 0.171746 kilome ters /day 2 
max 2 J 

In Figures 4 and 5 (pages 20 and 21) a (x, t) is plotted as a function of t 
for fixed values of x. 

Substituting (32a) and (32b) in the equations of motion for the satellite 
(5a) and (5b), we obtain: 


X - 2 cay - (1+2A) cu 2 x = K x x + 1 Kp [3 cos0t + 5 cos3rf,t] 

+ y K 3 [1 + 3 cos2(/)t] x (35a) 

y + 2a>x - (1-A) uu 2 y = £ K s [sintft + 5 sin30t] + j K 3 [sin2<Jt]x (35b) 

These are the linearized equations of motion for a satellite in the vicinity 
°f 1*2 j including both first and second order effects of solar perturbation. 



Table 1 


Perturbative Acceleration of Satellite at L P 


t (days) 


■ 


.00 

.17175 

-.00000 

.17175 

1.00 

.14915 

-.68500 x lO" 1 

.16413 

2.00 

.89765 x lO" 1 

-.11160 

.14322 

3.00 

.15478 x lO" 1 

-.11383 

.11488 

4.00 

-.46866 x lO" 1 

-.75687 x 10 _1 

.89022 x lO" 1 

5.00 

-.75944 x 10 _1 

-.13361 x 10- 1 

.77111 x lO" 1 

6.00 

-.64205 x 10- 1 

.47681 x lO -1 

.79973 x lO" 1 

7.00 

-.20666 x lO’ 1 

.82775 x 10" 1 

.85316 x lO" 1 

8.00 

.32836 x 10- 1 

.7780,1 x lO' 1 

.84446 x 10* 1 

9.00 

.70455 x 10- 1 

.34758 x 10- 1 

.78563 x 10- 1 

10.00 

.72738 x 10" 1 

-.29024 x 10- 1 

.78315 x lO" 1 

11.00 

.34490 x 10- 1 

-.87729 x lO' 1 

.94265 x lO -1 

12.00 

-.32981 x 10- 1 

-.11725 

.12180 

13.00 

-.10615 

-.10474 

.14913 

14.00 

-.15841 

-.53721 x 10 _1 

.16727 

15.00 

-.17044 

.17252 x lO -1 

.17131 

16.00 

-.13777 

.81807 x 10" 1 

.16023 

17.00 

-.72519 x 10- 1 

.11609 

.13688 

18.00 

.13288 x lO" 3 

.10807 

.10808 

19.00 

.57311 x lO -1 

.62245 x 10" 1 

.84611 x 10- 1 

20.00 

.76783 x lO" 1 

-.22258 x 10“ 3 

.76815 x lO -1 
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Table 2 


Perturbative Acceleration of Satellite Near L p 
(x = 1 kilometer) 




1 

BSflSSSflHH 

n 

.00 

-.33082 x 10- 1 

-.00000 

.33082 x 10' 1 

1.00 

-.55689 x 10- 1 

-.68683 x 10- 1 

.88423 x 10- 1 

2.00 

-.11511 

-.11193 

.16056 

3.00 

-.18946 

-.11426 

.22124 

4.00 

-.25186 

-.76126 x 10' 1 

.26311 

5.00 

-.28100 

-.13737 x 10- 1 

.28133 

6.00 

-.26930 

.47435 x 10- 1 

.27345 

7.00 

-.22579 

.82703 x lO' 1 

.24046 

8.00 

-.17228 

.77916 x lO' 1 

.18908 

9.00 

-.13463 

.35040 x 10' 1 

.13912 

10.00 

-.13230 

-.28627 x lO” 1 

.13536 

11.00 

-.17049 

-.87286 x lO’ 1 

.19154 

12.00 

-.23790 

-.11684 

.26505 

13.00 

-.31102 

-.10444 

.32809 

14.00 

-.36325 

-.53580 x lO" 1 

.36718 

15.00 

-.37527 

.17208 x 10" 1 

.37566 

16.00 

-.34262 

.81585 x lO" 1 

.35220 

17.00 

-.27741 

.11573 

.30058 

18.00 

-.20362 

.10764 

.23032 

19.00 

-.14770 

.61814 x lO" 1 

.16011 

20.00 

-.12828 

-.25762 x 10' s 

.12831 
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Table 3 


Perturbative Acceleration of a Satellite Near L ? 
(x = -1 kilometer) 


■■I 

1 

1 

1 1 

.00 

.37657 

.00000 

.37657 

1.00 

.35399 

-.68317 x 10" 1 

.36053 

2.00 

.29464 

-.11127 

.31495 

3.00 

.22041 

-.11341 

.24788 

4.00 

.15813 

-.75248 x 10- 1 

.17512 

5.00 

.12911 

-.12985 x lO" 1 

.12976 

6.00 

.14089 

.47927 x lO’ 1 

.14882 

7.00 

.18446 

.82847 x lO’ 1 

.20221 

8.00 

.23795 

.77686 x lO" 1 

.25031 

9.00 

.27555 

.34477 x 10" 1 

.27769 

10.00 

.27778 

-.29422 x 10" 1 

.27933 

11.00 

.23947 

-.88172 x 10- 1 

.25519 

12.00 

.17194 

-.11766 

.20834 

13.00 

.98716 x 10- 1 

-.10505 

.14415 

14.00 

.46425 x 10- 1 

-.53863 x 10- 1 

.71109 x 10 _1 

15.00 

.34389 x 10" 1 

.17297 x 10" 1 

.38494 x 10 _i 

16.00 

.67081 x 10" 1 

.82030 x 10- 1 

.10597 

17.00 

.13237 

.11645 

.17630 

18.00 

| 

.20628 

.10851 

.23307 

19.00 

.26232 

.62676 x lO" 1 

.26970 

20.00 

1 

.28185 

-.18754 x lO' 2 

.28186 
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Table 4 


Perturbative Acceleration of a Satellite Near L ? 
(x = 0.5 kilometers) 


t (days) 

a (kilome- 
ters /day 2 ) 

a (kilome- 
t?rs/day 2 ) 

a (kilome- 
ters/day 2 ) 

.00 

.69322 x lO" 1 

-.00000 

.69332 x lO" 1 

1.00 

.46732 x lO" 1 

-.68591 x 10" 1 

.82998 x 10" 1 

2.00 

-.12675 x 10" 1 

-.11177 

.11248 

3.00 

-.86988 x 10- 1 

-.11404 

.14343 

4.00 

-.14936 

-.75907 x 10" 1 

.16755 

5.00 

-.17847 

-.13549 x 10- 1 

.17898 

6.00 

-.16675 

.47558 x 10" 1 

.17340 

7.00 

-.12323 

.82739 x lO’ 1 

.14843 

8.00 

-.69724 x 10- 1 

.77859 x lO -1 

.10451 

9.00 

-.32090 x lO" 1 

.34899 x 10" 1 

.47410 x 10 _1 

10.00 

-.29782 x lO" 1 

-.28825 x lO" 1 

.41447 x lO" 1 

11.00 

-.68000 x 10" 1 

-.87508 x 10 _1 

.11082 

12.00 

-.13544 

-.11705 

.17901 

13.00 

-.20859 

-.10459 

.23334 

14.00 

-.26083 

-.53651 x 10" 1 

.26629 

15.00 

-.27285 

.17230 x lO" 1 

.27340 

16.00 

-.24019 

.81696 x 10 _1 

.25371 

17.00 

-.17496 

.11591 

.20987 

18.00 

-.10114 

.10785 

.14786 

19.00 

-.45194 x 10 _1 

.62030 x 10- 1 

.76748 x 10 _1 

20.00 

-.25750 x lO -1 

-.24010 x 10" 2 

.25862 x 10- 1 
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Table 5 


Perturbative Acceleration of a Satellite Near L s 
(x = -0.5 kilometers) 


jHH 

■ - ■ 1 

HSPhHHI 

mMMBmBmm 

.00 

.27416 

.00000 

.27416 

1.00 

.25157 

-.68409 x 10- 1 

.26071 

2.00 

.19220 

-.11143 

.22217 

3.00 

.11794 

-.11362 

.16377 

4.00 

.55632 x 10- 1 

-.75468 x 10- 1 

.93756 x lO -1 

5.00 

.26583 x 10- 1 

-.13173 x lO -1 

.29668 x 10 _1 

6.00 

.38344 x 10- 1 

.47804 x 10" 1 

.61282 x lO -1 

7.00 

.81895 x 10- 1 

.82811 x 10- 1 

.11647 

8.00 

.13539 

.77743 x 10- 1 

.15613 

9.00 

.17300 

.34617 x lO -1 

.17643 

10.00 

.17526 

-.29223 x 10- 1 

.17768 

11.00 

.13698 

-.87950 x lO" 1 

.16279 

12.00 

.69479 x 10- 1 

-.11746 

.13647 

13.00 

-.37179 x 10~ 3 

-.10490 

.10496 

14.00 

-.55993 x lO" 1 

-.53792 x lO" 1 

.77645 x lO" 1 

15.00 

-.68025 x lO" 1 

.17275 x 10- 1 

.70184 x 10- 1 

16.00 

-.35343 x 10" 1 

.81919 x 10- 1 

.89218 x 10- 1 

17.00 

.29926 x 10 _1 

.11627 

.12006 

18.00 

.10380 

.10829 

.15000 

19.00 

.15982 

.62461 x 10- 1 

.17159 

20.00 

.17932 

-.20506 x lO" 3 

.17933 
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Table 6 


Perturbative Acceleration of a Satellite Near L^, 
(x = 50 kilometers) 


1 


■ 

^ l 

a (kilome- 
ters /day 2 ) 

.00 

-.10070 x 10 2 

-.00000 

.10070 x 10 2 

1.00 

-.10093 x 10 2 

-.77643 x 10' 1 

.10093 x 10 2 

2.00 

-.10154 x 10 2 

-.12825 

.10155 x 10 2 

3.00 

-.10231 x 10 2 

-.13503 

.10232 x 10 2 

4.00 

-.10297 x 10 2 

-.97640 x 10- 1 

.10297 x 10 2 

5.00 

-.10329 x 10 2 

-.32156 x 10 _1 

.10329 x 10 2 

6.00 

-.10319 x 10 2 

.35398 x lO -1 

.10319 x 10 2 

7.00 

-.10277 x 10 2 

.79195 x 10" 1 

.10277 x 10 2 

8.00 

-.10223 x 10 2 

.83562 x 10- 1 

.10223 x 10 2 

9.00 

-.10184 x 10 2 

.48833 x 10- 1 

.10184 x 10 2 

10.00 

-.10179 x 10 2 

-.91466 x 10- 2 

.10179 x 10 2 

11.00 

-.10215 x 10 2 

-.65596 x 10- 1 

.10215 x 10 2 

12.00 

-.10279 x 10 2 

-.96813 x 10- 1 

.10279 x 10 2 

13.00 

-.10350 x 10 2 

-.89644 x 10" 1 

.10350 x 10 2 

14.00 

-.10400 x 10 2 

-.46655 x 10- 1 

.10400 x 10 2 

15.00 

-.10412 x 10 2 

.15024 x lO" 1 

.10412 x 10 2 

16.00 

-.10380 x 10 3 

.70681 x 10- 1 

.10380 x 10 2 

17.00 

-.10317 x 10 2 

.98049 x 10" 1 

.10317 x 10 2 

18.00 

-.10246 x 10 2 

.86339 x lO" 1 

.10246 x 10 2 

19.00 

-.10193 x 10 2 

.40697 x 10- 1 

.10193 x 10 2 

20.00 

-.10177 x 10 2 

-.19746 x 10- 1 

.10177 x 10 2 


- 18 - 




Table 7 


Perturbative Acceleration of a Satellite Near L 3 
(x = -50 kilometers) 


■■1 

H 


a (kilome- 
ters /day 2 ) 

.00 

.10413 x 10 S 

.00000 

.10413 x 10 2 

1.00 

.10391 x 10 2 

-.59358 x 10- 1 

.10391 x 10 2 

2.00 

.10334 x 10 2 

-.94946 x 10" 1 

.10334 x 10 2 

3.00 

.10262 x 10 2 

-.92637 x 10" 1 

.10263 x 10 2 

4.00 

.10203 x 10 s 

-.53735 x lO" 1 

.10203 x 10 2 

5.00 

.10177 x 10 2 

.54331 x 10- 2 

.10177 x 10 s 

6.00 

.10191 x 10 2 

.59964 x 10" 1 

.10191 x 10 2 

7.00 

.10235 x 10 2 

.86355 x 10’ 1 

.10236 x 10 2 

8.00 

.10289 x 10 2 

.72040 x lO" 1 

.10289 x 10 2 

9.00 

.10325 x 10 2 

.20683 x 10 _1 

.10325 x 10 2 

10.00 

.10325 x 10 2 

-.48902 x lO" 1 

.10325 x 10 2 

11.00 

.10284 x 10 2 

-.10986 

.10284 x 10 2 

12.00 

.10213 x 10 2 

-.13769 

.10214 x 10 2 

13.00 

.10137 x 10 2 

-.11985 

.10138 x 10 2 

14.00 

.10083 x 10 2 

-.60787 x lO -1 

.10084 x 10 2 

15.00 I 

.10071 x 10 2 

.19481 x 10 _1 

.10071 x 10 2 

16.00 

.10105 x 10 2 

.92933 x 10' 1 

.10105 x 10 2 

17.00 

.10172 x 10 2 

.13412 

.10173 x 10 2 

18.00 

.10249 x 10 2 

.12980 

.10250 x 10 2 

19.00 

.10308 x 10 2 

.83794 x 10- 1 

.10308 x 10 s 

20.00 

.10330 x 10 2 

.15295 x 10" 1 

.10330 x 10 2 
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Figure 5. Perturbative Acceleration of Satellite x20 



3. ANALYTICAL SOLUTION OF THE EQUATIONS OF MOTION 

Since no closed form solution is known to exist for equations (35a) and 
(35b), we will develop approximate analytic solutions of first and second 
order. The constant K 3 , which appears on the right hand side of eqations 
(35a) and (35b), is far smaller than K x and Kg as can be seen from the numerical 
values of these constants given in this report. Therefore, it is logical to 
express both x and y in powers of K 3 and then solve the equations of motion by 
the method of regular perturbations. 

For conformity, let e = K 3 . Then: 

x = x x + exg + e 2 x 3 + ... (36a) 

y = Yi + ey 2 + e 2 y3 + • • • (36b) 


where the subscripts denote the order of solution. 


The initial conditions will be taken as 


x i (0) 

x x (0) 

yi(0) 

y, (0) 

* 3 ( 0 ) 

*s(0) 

y 2 (°) 

y 8 (0) 


x(0) 


x(0) 


y(0) 


y(0) 


x 3 (0) = ., 

.. = 0 

X • 

CO 

✓**N 

o 

V 

II 

• 

.. = 0 

y 3 (0) = ., 

,. = 0 

y 3 (0) ■ •• 

.. = 0 


(37) 


In this section we will derive the solutions corresponding to €° and e 1 . There 
is little point in going to higher order terms since the equations of motion are 
linear lized, and are valid only in the vicinity of L a . 
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Substituting the expansions (36a, b) in (35a, b) we obtain two sets of differen- 
tial equations corresponding respectively to e° and e 1 . 

For e° : 

• 1 

xj - 2y<yi - Kg [3 cos 0 t + 5 cos30t] (38a) 

y x + 20 ^ - 8 y x = ^ Kj, [sin 0 t + 5 sin3</>t3 (38b) 

where or = (l+ 2 A){t ? 2 + K z = 0.184100 

and 8 = (1-A)f!> 3 »-0. 115282 

Similarly, for e 1 : 

x g - 2(i.)y s - axg = ~ [1 + 3 cos20tlx 1 (39a) 

y p + 2'jtXg - 8 y 3 = j [sin20t]x 1 (39b) 


The complete second-order solution will then be given by 

x = x x + K 3 Xj, (40a) 

y = yi + K 3 y s (40b) 


3.1 FIRST ORDER SOLUTION (e°) 

This section will deal strictly with the system (38a, b) . We first solve 
the set of homogeneous equations, which gives rise to a fourth-order characteris- 
tic equation: 

D 4 + D 2 (4f' 2 - cv - 3) + aP = 0 (41) 

which has two equal and opposite real roots denoted by + p , and two equal and 
opposite imaginary roots denoted by + iQ. 


P - [- (2-A )(! 3 + K, + /^‘ (9A 2 -8A) + ( 6 A- 8 ) + kJ" 1 ^ 


= 0.301427 


Q = y= r( 2 -A), 3 - K x + /u 4 (9A 3 -8A) '+"'u?k 1 ( 6 A- 8 ) + K , 2 ] 


./ 3 


= 0.483308 


(42) 


(43) 
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The homogeneous solution is, therefore: 


■ H 


- H 


= A 1 sinQt 

+ Ap 

cosQt + A 3 

sinh 

pt + a 4 

cosh 

pt 

(44a) 

= B x sinQt 

+ B P 

cospt + B s 

sinh 

pt + b 4 

cosh 

pt 

(44b) 


The B's are related to the A's as follows: 

B x ~ ~v Ap 
B 2 = y A x 
B a =5A 4 

B,i = 8 A 3 

where : 


(45) 


Y = 2 ^ rn 3 + (l+2A).n 3 + Kj 1 = 1.878986 

and : 

(46) 

& = r-t- r p 3 - ( l+2A)u) 3 + K,1 = -3 . 631650 

2 cup 1 - 

(47) 

The particular solution is found to be : 


X X (t) - A 5 COSQt + Ag COS 3rf t 

(48a) 

Yi (t) = B 5 sin<At + sin3r^t 

(48b) 

where : 


K? 2'jjri - 3 (^ 2 +8) 

5 4 (o ,+/ * 2 ) (P+P 3 ) - 4 u) 2 ? 3 

(49a) 

5K S 6(juP - (9- i2 +3) 

e 4 (or^Q 3 ) (|3 + 9* 2 ) - 36 n s 0 2 

(49b) 

- - 2^ -K + 

(49c) 

B « - - isr|K a + (0+9^)*,] 

(49d) 


The remaining four constants can be determined in terms of the four given 
initial conditions. Namely: 


A = Py.CQ) - Sx(0) 

1 py - 05 




- &P [x(0) - A, - A,J + $ [B s + 3 B 6 ]- y(Q) 


f^Y 


+ p6 


(50a) 

(50b) 
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(50c) 


A3 ~ p [ X (^) “ OA\ 1 

A 4 = x(0) “ Ap - A 5 - A 6 (50d) 

The general first order solution is given by: 

x 1 (t) = A l sinQt + Ag cosQt + A 3 sinh pt + A^ cosh pt + A 5 cos0t 

+ Ag cos30t (51a) 


y 1 (t) = B x sinQt + B 2 cosQt + b 3 sinh pt + B 4 cosh pt + B 5 sin^t 

+ B 6 sin30t (51b) 

where all the constants have been defined either explicitly or in terms of the 
initial conditions. It is evident that this solution is unbounded due to the 
presence of the hyperbolic functions. Even for small values of t the exponential 
terms are larger than the sinuosoidal terms. Consequently, this solution 
demonstrates little, if any, periodic behavior even in the initial phase of the 
trajectory. The source of the instability can be eliminated, however, by a proper 
choice of initial conditions. That is, we want the coefficients of the hyperbolic 
terms A 3 , A 4 , B 3 , B 4 , to vanish. As a matter of fact, we must only require that 
A 3 = A 4 = 0 since B 3 and B 4 are multiples of A* and A 3 respectively (see (45)). 

By letting A 3 = 0 and A 4 = 0 (in (50c) and (50d)) we determine new expressions 
for A x and A^ , namely: 



(52a) 


Ap = x(0) - A 5 - A 6 (52b) 

Substituting the above in (50a) and (50b) we come up with the necessary relation- 
ship among the initial conditions: 


x(0) = 2 y(0) (53a) 

y(0) = 0(B s +3B e ) + yQ (A 5 +A 6 ) - yQx(O) (53b) 

These two relations can be implemented in practice during the injection into an 
orbit around Lg . Thus we have a periodic first order solution, namely: 

x i (t) = Ai sinQt + Ap cosQt + A 5 costft + A s cos30t (54a) 
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y x (t) = B x sinQt + B s cosQt + B 5 sin0t + B 6 sin30t 


(54b) 


where : 

Ab = -.09389 
Ag = -.72532 
B 5 = 0.43798 
Bg = -1.09595 

and: 

A l = y(0)/y = 0.53220 y(0) 

Ag = x( 0)-A 5 -A 6 = 0.81922 + x(0) 

Bi = -yA 3 = -1.53930 - 1.87896 x(0) 

Bg = y(0) 

Using the first order solution we can then compute, in addition to the trajec- 
tory, the velocity, acceleration, range to L s , range rate, and flight path angle. 
Namely: 


v(t) = /x x 3 (t) + y x 2 (t) 

(55) 

a(t) = A a (t) + yi 2 (t) 

(56) 


(57) 

(58) 

(59) 

3.2 SECOND ORDER SOLUTION (e 1 ) 

Having found the first-order solution, we can proceed to determine the 
second-order solution using equations (39a, b) . Substituting x x (t) from equa- 
tion (54a) into equations (39a, b) , we obtain the following system: 

Xg - 2(oy s - q qc, = j [ 1+3 cos20tl[A x sinQt + Aj, cosQt + Ag cos0t 

+ A g cos30t] (60a) 
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9 2 

y 2 + 2 ojXg - 3y s = — [sin20t][A 1 sinQt + cosQt + A 5 cos0t 

+ A 6 cos30t] (60b) 

Expanding the right hand sides of the above equations and separating all 
cross products of trigonometric functions, we obtain a simplified version which 
is easily solved. 

* 1 1 1 

Xg - 2(ry p - OXg = — A x sinQt + — A s cosQt + [5 A s + 3 A 6 ] cos0t 

+ ^ [2 Ag + 3 Ajl cos30t + ^ A s cos50t 
+ £ A x sin (20+Q) t - ^ A x sin (20-Q)t 

+ ^ Ap cos (20+Q)t + ^ AgCOS (20-Q)t (61a) 

• 3 3 3 

y ? + 2a>Xg - gy s = ^ [Ag-Ae] sin0t + ^ A 5 sin3^>t + ^ Ae sin50t 

3 3 

+ £ Aj cos (20-Q)t - A x cos (20+Q)t 

+ ^ Ag sin ( 20+Q) t + ^ A 2 sin (20-Q)t (61b) 

The homogeneous solution is unchanged from the first-order solution except for 
the coefficients. That is: 

Xg H (t) = Aj sinQt + Ag' cosQt + A 3 ' cosh pt + A*' sinh pt (62a) 

y-5„(t) = Bi' sinQt + B s / cosQt + B3 cosh pt + B 4 sinh pt (62b) 

H. 

where ; 

Bi = -yAg 

= VAi (63) 

B3 - SA^ 

B 4 ' = 5Ag 
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The particular solution is found to be: 


Xj,p(t) = Cj COS0t + C; 

+ C 5 sin (20' 
+ C e sinQt + 

y 8 (t) = Di sin 0 t + D. 

+ D s sin (20 
+ D s sinQt + 


cos30t + C 3 cos50t + C 4 
•Q)t + C 6 cos (20+Q)t + C- 
C g cosQt 

sin30t + D a sin50t + D 4 
•Q)t + D s cos ( 2 0 +Q) t + D. 
D 9 cosQt 


sin (2^+Q)t 
cos ( 20 -Q)t 

(64a) 

sin (2rA+Q)t 
cos ( 20 -Q)t 

(64b) 


It should be noted that the sinuosoidal part of the homogeneous solution 
appears in the particular solution without causing resonance. This is due to 
the particular way in which the differential equations are coupled. The con- 
stants are defined below. 


3/2 a* (As-Afi) - 1/4 (0 3 +B)(5A 5 + 3A«) 
(0 2+ a?) (* 3+ B) - 4(j^ 0 3 

(65a) 

9/2 uxZ> A k - 1/4 (9<* 3+ B)(3A K + 2As) 
(9* 3 +a) (90 3 +g) - 36a; 3 0 1 3 

(65b) 

3/4 Ag (lOa>0 -250 s -B) 
(25^+00 (250 3 +B) - lOOto^ 3 

(65c) 

3/4 A!r(2rt+Q) 3 -2co (2<ft+fi) + Bl 
f ( 2 *+Q) 3 +0' 1 [ ( 2 *+Q) 3 +8 1 - 4<* 3 (2 ,-a+Q) 8 

(65d) 

3/4 Air( 2 ^-Q) 3 -2aj (20-Q) + B 1 

r(20+Q) 3 +Q'] r (2'*-Q)‘ i +B^-4<.>' 3 (2^-Q) 3 

(65e) 

3/4 A?f(2rt+Q) 2 - 28) (2*+fi) + Bl 
4 2 v+O)^ - f ( 2 0+Q) 2 +0 ] [ ( 2irt+Q) 3 +B ] 

( 65 f ) 

3/4 A~r(2tf-Q) 3 - 2 co ( 2 +-Q) + Bl 
4u) 3 (20-o)s-[(2^-n) s +Q'l[(2M-Q)2+B] 

(65g) 

1 (8+Q 2 ) Ai 

2 [(a+Q s )(B+n 3 )- 4 

(65h) 
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c e = . £ 


(g+o 2 )^ 


Di = 


D„ = 


D, = 


D 4 = 


D* = 


2 [(a+O^Kg+O 2 ) + 4 u£Q3] 

2(r0Ci + 3/4 (As-Afi) 

(<f + 3) 

6qj^C? + 3/4 As 
(90 s + 3) 

lOg)0C3 + 3/4 Ar 
25<f + g 


D 7 = 
D 8 = 
D 9 = 

A/ = 

K = 


(2(f) +0) 2 

+ 3 

2u/20-O) C 7 - 

f 3/4 A? 

(20-0)3 

+ 3 

2fo (20+0) C 4 

+ 3/4 Ai 

(20 + 0) 2 

+ 3 

2u) (20-0) Cs 

- 3/4 A) 

(2^ - 0) 2 

+ 3 

1/2 As + (a+O 2 ) C 9 


2cuQ 

1/2 Ai + (cy+Q 2 ) Cb 

2ojQ 


|5[C 4 (20+0) + C 5 (20-fi) + C S Q] - p [D 6 + D 7 + D 9 ]| 
{d x 0 + 3 D ? 0 + 5 D 3 0 + (20+0) D 4 + (20-0) D 5 + D S Q 


- Sp [C x + C s + C 3 + C 6 + C 7 + C 9 ]j- 

A s = - [C x + c 2 + C 3 + C p + C 7 + C e ] - Ag 

A 4 = - | [D s + D 7 + D 9 ] - | A x 
The second order solution is thus given by: 

x 2 ( fc ) = ( A i + c e) sin 0t + (Aj, + C 9 ) cosQt + Ag cosh pt + A 4 ' sinh pt 

+ C x cos0t + C 2 cos30t + C 3 cos50t + C 4 sin (2^+0) t 

+ C s sin (20-Q) t + C s cos (20+0) t + C 7 cos (20-Q)t 


(65i) 

(66a) 

(66b) 

(66c) 

(66d) 

(66e) 

(66f) 

(66g) 

(66h) 
( 6 6i ) 

(67a) 


(67b) 

(67c) 

(67d) 


(68a) 
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y p (t) = (b i + D s ) sinQt + (bJ + D 9 ) cosQt + B 3 cosh pt + B 4 sinh pt 


+ D x sinQt + D p sin3^t + D 3 sin50t + D 4 sin (20+Q)t 

+ D 5 sin (20-Q)t + D 6 cos (20+Q)t + D ? cos (20-Q)t (68b) 

In the case of this solution the hyperbolic terms are left intact to 
demonstrate the increasing instability of the solution. 

3.3 COMPLETE SECOND ORDER SOLUTION 

The complete second order solution is given by: 

x(t) = x x (t) + K 3 x p (t) (69a) 

y(t) = y x (t) + K 3 y P (t) (69b) 

where x 1 (t), x p (t), y 1 (t) and y ? (t) are given by (54a), (54b), (68a) and (68b), 
respectively. 

Thus the complete solution becomes: 

x ( t ) = [A x + K 3 (Ai+Cg)] sinQt + [Ag + K s (Ag'+Cg)] cosQt 

+ K 3 A 3 cosh pt + K3A4 sinh pt + [Ag+KsCi] cos0t + [A 6 -4C 3 C ? ] cos30t 
+ K 3 C 3 cos50t + K 3 C 4 sin (20+Q)t + K 3 C B sin (2Q-Q)t 

+ K 3 C S cos (2*+Q)t + K 3 C 7 cos (20-Q)t (70a) 

y(t) = [B 1 + K 3 (Bj'+Dg)] sinQt + [B a -tK 3 (Bp'+Dp)] cosQt 

+ K 3 B 3 cosh pt + K 3 B 4 sinh pt + [Bg+KsD]^ ] sin0t + [B P +K 3 D ? ] sin30t 
+ K 3 D 3 sin5^t + K 3 D 4 sin (20+Q)t + K 3 D 5 sin (20-Q)t 

+ K 3 D p cos (20+Q) t + K 3 D7 cos (20-Q)t (70b) 

The numerical values of the coefficients will be given here for a typical 
trajectory originating at L s . 
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Since the initial conditions chosen are periodic only as far as the first order 
solution is concerned, we will rename them quasi-periodic for the complete 
solution. Thus, the quasi-periodic initial conditions for a start at L P are: 

x(0) = 0 
y(0) = 0 


x(0) 

= 

0 







7(0) 

= 

-.13750 







Then : 









Ai 

= 

0 

K 

B 

0 

B2 

= 

1.9471 

As 

= 

0.81922 

B3 

= 

0 

d 3 

B 

0.89387 

A* 

s= 

-0.09389 

B 4 

b 

18.2830 

d 4 

= 

4.5579 

Ag 

= 

-0.72532 

Ci 

= 

3.6139 

d 5 

= 

-1.6151 

B i 

B 

-1.53929 


= 

1.6983 

Ds 

= 

0 

b 3 

- 

0 

c 3 

= 

0.74566 

d 7 

B 

0 

B5 

= 

0.43798 

c 4 

B 

0 

Da 

= 

-0.92132 

Be 

= 

-1.09595 

Cg 

= 

0 

d 9 

= 

0 

Ai 

= 

0 

c 6 

S 

-3.9235 

k 3 

= 

0.29525 11 

Ap 

=3 

4.6669 

c 7 

= 

-1.2766 

Q 

= 

0.4833081 

A3 

= 

-5.0343 

Ca 

= 

0 

0 

= 

-0.2128 

A4 

= 

0 

C 9 

s 

-0.49033 




Bi' 

b 

-8.7691 

Di 

= 

1.7124 





3.4 TABULATIONS AND GRAPHS OF SAMPLE TRAJECTORIES 

The results of our computations will be tabulated and plotted in the follow- 
ing pages. First, the first order solution will be presented with periodic ini- 
tial conditions for a few different cases. Then the complete second order solu- 
tion, with quasi-periodic initial conditions, will be given. The instability of 
the motion will be demonstrated in this solution. However, it may be noted, the 
divergence of the motion commences only after 23 days for a motion initiating 
at Lg , with quasi-periodic initial conditions. Thus the rate of growth of the 
divergence from the periodic motion is rather slow and probably could be 
corrected with the aid of an onboard engine. 
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First Order Solution with Periodic Initial Conditions 
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Table 8 (Continued) 

First Order Solution with Periodic Initial Conditions 


I 
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Y (km) 



Figure 6. Trajectory of Satellite Around l _2 
First Order Solution 


Periodic 

Initial Conditions 
x = 1.61 km 
y — 1 .61 km 
x = 0.414 km/ day 
y = -1 .60 km/ day 
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ure 7. Velocity of Satellite Versus 
First Order Solution 





Figure 8. Acceleration c 
First Ordei 
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Figure 9. Range of Satell 
First Order 






First Order Solution with Periodic Initial Conditions 
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Table 9 (Continued) 

First Order Solution with Periodic Initial Conditions 
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Figure 11. Velocity of Satellite Versus Time 
First Order Solution 










Figure 14, Range Rate of Satellite Versus Time 
First Order Solution 
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Table 10 (Continued) 

First Order Solution with Periodic Initial Conditions 
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Figure 22. Velocity of Satellite Versus Time 
Complete Second Order Solution 






Figure 23. Acceleration of Satellite Versus Time 
Complete Second Order Solution 
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Figure 24. Range of Satellite from L 2 Versus Time 
Complete Second Order Solution 




= 0.8 km 
= 0.209 km/ day 
= -0.869 km/ day 
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